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A.  STATEMENT  OF  THE  PROBLEM  STUDIED 


1.  Motivation 


The  importance  of  conservation  laws  is  widely  recognized  in  many  fields  of  physics 
and  mechanics.  Conserved  quantities  such  as  ma.ss.  energy,  momentum  and  angular  mo¬ 
mentum  are  fundamental  in  classical  and  modern  physics.  They  provide  insights  into  the 
behavior  of  various  systems,  ranging  from  dynamic  of  rigid  bodies  to  continuum  physics. 
In  the  theory  of  elasticity,  conservation  laws  represented  as  path-independent  integrals 
(such  as  J,  L  and  M)  have  captured  the  interests  of  many  researchers  (see.  for  example, 
[7),  [S]  ,[9|  ,[10]  and  [14]).  These  integrals  contribute  significantly  to  the  study  of  fracture 
and  defects.  Conservation  laws  may  also  be  useful  in  the  improvement  of  algorithms  in 
numerical  solutions,  the  establishment  of  global  ('xistence  and  uniqueness  theorems,  the 
study  of  stability  of  solutions,  the  investigation  of  pha.se  transformation  and  the  reduc¬ 
tion  of  the  order  of  the  governing  equations  for  different  systems  considered.  Since  there 
are  numerous  applications  for  conservation  laws,  it  is  of  considerable  interest  to  develop 
methodologies  for  obtaining  these  laws  and  to  apply  them  to  \ariou.s  sj'stems  of  practical 
interest,  including  dissipative  systems. 


2.  General  Background 


Mathematically,  a  conservation  law  of  a  physical  system  with  four  independent  vari¬ 
ables,  X,  y.  z  and  f,  for  example,  is  an  (’(|uatiou  of  the  form 

Div  P  =.  D.  P,  D::P„  +  D,P,  +  D,P,  ( 1 ) 

where  P  =  {P^,  Py.  P..  Pi)  is  a  vector  function  that  depends  on  tlu'  independent  varialdes. 
the  dependent  variables  and  dcrivariv<'s  of  the  cL'pcmdeut  varial^les  of  th('  s\'stem.  Here 
D,  designryes  a  total  derivative  operator.  Physically,  a  conservation  law  states  that  the 


rate  of  change  of  Pt  inside  any  spatial  domain  is  ec[ual  to  the  not  Hux  of  the  “currents” 
{Pxi  Pyi  Pz)  through  the  surface  of  the  domain.  In  the  special  ciise  of  systems  with  only 
one  independent  variable,  conservation  laws  lead  to  qm'uitities  which  are  constant  on  the 
solution  manifold. 

Based  on  variational  principles,  Noether's  first  theorem  provides  a  method  to  estab¬ 
lish  conservation  laws  for  Lagrangian  systems,  i.e..  for  systems  possessing  a  Lagrangian 
function  and,  therefore,  governed  by  the  associated  Euler-Lagrange  equations.  While  pro¬ 
viding  a  condition  for  e.xistence  of  conservation  laws.  Noether’s  theorem  by  itself  does  not 
offer  an  explicit  method  for  constructing  divergence-free  expressions.  However,  if  one  com¬ 
bines  Noether’s  theorem  with  Lie  group  theory,  a  systematic  procedure  for  constructing 
conservation  laws  for  Lagrangian  systems  becomes  available  [2j,  [15]  and  [16]. 

While  providing  a  direct  procedure  for  obtaining  conservation  laws,  Noether’s  theo¬ 
rem  is  applicable  only  to  Lagrangian  systems.  In  1921,  Bessel-Hagen  extended  Noether’s 
theorem  by  introducing  the  concept  of  divergence  symmetries,  however,  his  generalization 
also  operated  only  in  realm  of  Lagrangian  systems 

Several  isolated  attempts  have  previously  Ix-xm  made  with  the  aim  of  deriving  consci- 
vation  laws  for  dissipation  systems.  Most  of  these  start  Iw  trying  to  reduce  the  problem 
to  a  variational  framework  in  which  Noether’s  theorem  is  applicable.  For  example,  Jiang 
[12]  and  [13]  has  developed  conservation  laws  in  viscoelastostatics  and  viscoekrstodynamics 
by  employing  first  Laplace’s  transform  and  then  applying  .Noether's  first  theorem  to  the 
transformed  eciuations.  The  conservation  laws  thus  ol)tained  involve  convolution  and  when 
expressed  in  the  time-domain  they  therefore  depend  on  the  total  history  of  the  material 
evolution.  Later,  Caviglia  and  Morro  [3]  applieil  the  same  procedure  with  the  Bessel-Hagen 
extension  of  Noetlier  s  tlieorem  to  arrive  at  more  general  results.  Soon  afterwards  [4].  [5], 
they  reexamined  their  results  in  view  of  a  differi'iit  procedure,  called  direct  approach, 
which  has  been  developed  b\'  .Nrens  [l]  in  connection  wirh  Maxwell  equations.  The  direi  t 
approach  consists  simply  m  integrating  eijuarion  i  1 )  ^uiijecr  to  the  constraint  provided  b}' 
the  governing  e([uation  of  the  systc-m. 

In  the  rCcalm  oi  C(?nser\'ation  laws  ior  dis'^mativ'e  s\'stenis.  u  is  worth  mi'iitii'iuing  the 


ficient  condition  for  the  derivation  of  the  conservation  laws. 


In  a  recent  work  [11],  the  present  authors  have  i)roposed  a  novel  approach,  called 
the  Neutral  Action  method,  to  derive,  in  a  highly  systematic  fashion,  conservation  laws  for 
dissipative  systems  governed  by  ordinary  as  well  as  partial  differential  ecpiatioas.  Moreover, 
the  conditions  provided  for  the  existence  of  con.sorvation  laws  are  necessary  as  well  as 
sufficient. 


3.  Statement  of  the  problem 


The  problem  consists  in  applying  the  Neutral  Action  method  in  order  to  derive  con¬ 
servation  laws  for  viscoelastic  bodies.  In  contrast  with  the  previous  attempts,  emphasis 
will  be  put  in  constructing  conservation  laws  in  the  space-time  domain  and  the  viscoelastic 
models  considered  are  therefore  given  in  terms  of  rate  e([uations. 


B.  SUMMARY  OF  THE  RESULTS  ACHIE\'ED 


The  “Neutral  Action”  method  is  ajjplicnble  to  dissipative  as  well  as  nondissipative 
systems. 

specialized  to  Lagransiaii  systems,  it  reduces  to  Noether’s  first  theorem  as  ex¬ 
tended  by  Bessel-Hagen.  By  applying  the  “Neutral  Action”  method  we  were  Liierefore  able 
to  obtain  in  a  unified  and  systematic  hishittn  all  th<'  aforementioned  resrdts.  In  addition, 
conservation  laws  in  the  space-time  domain  were  derived  for  a  one-dimensional  viscoelas¬ 
tic  bar.  The  models  considered  first  were:  Voigt,  Kt-lvin  and  Maxwell  elements.  Then  we 
considered  the  more  general  case  where  the  constitutive  equation  is  given  by 

Q{d,)a  =  R{d,)e  (2) 

where  Q  and  R  are  polynomials  in  the  operator  di.  and  a  and  e  designate  the  stress  and 
strain,  respectively.  For  the  quasistatic  case,  we  w('re  able  to  show  that  the  conserN'ation 
laws  are  given  by 

■Fr  =  )  -  C  -  .ij  (3) 

Pi  ~  <AOx)  (4) 

where  6  is  defined  by  0  =  R(dt)u,  u  being  the  di.splacf'ini'nt  field  and  g  and  h  are  arbitrary 
functions,  while  C  is  an  arbitrary  constant.  .A  prime  designates  diffcrt'ntiatiou  with  res]:)ect 
to  the  arguments. 

By  suitable  choices  of  g  and  /i,  different  conservation  laws  can  be  obtained.  One  of 
them  is  the  statement  of  energy  coirservatiou. 

A  generaliziition  to  two-dimension;tl  linear  i^.otropie  viscoelastic  mati'rials  has  been 
carried  out  for  the  Voigt  model. 

One  of  the  conservation  laws  obtained  is  given  by 

Px  -  +  iirr,,-,,  +  i] 

P'/  “  4-  I  -f  >'0,1,1,  -r  no ,  ,j 


P*  —  —A  {2Vy'  +  (Ct  +  2/3)(ti  j;U  j;/  +  +  Ol{v^yU  :^t  + 


+  0{U,yU  ,yt  ^  ,z^  ,yt  ^,y^.zt)\ 

where  A  is  an  arbitrary  constant,  W‘  is  the  elastic  strain  energy  and  c 
cients  of  dissipation.  One  can  verify  that  this  conservation  law  is  nothine 
of  energy  conservation. 


V,  3  are  the  coeffi- 
;  but  the  statement 
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